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Abstract 


The National Aeronautics and Space Administration 
(NASA) is conducting research with the goal of 
enabling safe improvements in the capacity of the 
nation’s air transportation system. The wake- vortex 
upset hazard is an important factor in establishing the 
minimum safe spacing between aircraft during landing 
and take-off operations, thus impacting airport capacity. 
A batch simulation study was conducted to assess the 
sensitivity of various safe landing criteria in the 
development of an acceptable wake encounter 
boundary. A baseline six-degree-of-freedom simulation 
of a B737-100 airplane was modified to include a wake 
model and the vortex-induced forces and moments. The 
guidance and control input for the airplane was 
provided by an auto-land system. The wake strength 
and encounter geometry were varied. A sensitivity 
study was also conducted to assess the effects of 
encounter modeling methods and accuracy. 

Nomenclature 

C { rolling moment coefficient 

rolling moment due to sideslip coefficient 

C L lift coefficient 

lift coefficient change with angle of attack 

yawing moment due to sideslip coefficient 

C Y sideforce coefficient 

C y ^ sideforce due to sideslip coefficient 

h„ altitude of the wake above the runway, ft 
r c vortex core radius, ft 

r radius from vortex center, ft 
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v velocity component along Y axis, ft/s 



vortex tangential velocity, ft/s 

w 

velocity component along Z axis, ft/s 

X 

longitudinal axis 

y 

coordinate along Y axis, ft 

y* 

lateral coordinate of the wake, ft 

Y 

lateral axis 

z 

coordinate along Z axis, ft 

Z 

vertical axis 

Sw 

roll control power ratio 

r 

vortex circulation strength, ft 2 /s 

Subscripts 

auto 

referenced to the autopilot system 

b 

reference to body axis system 

L 

reference to left or port vortex 

max wheel with maximum wheel input 

no wake 

with no wake present 

P 

reference to point P 

rwy 

reference to runway axis system 

R 

reference to right or starboard vortex 

w 

wake-induced component 

wake 

with wake present 


Introduction 


Many of today’s major airports are capacity limited, 
leading to increased airport congestion and delays. The 
ability to relieve the congestion through airport 
expansion or new airport construction is limited and 
increasingly difficult. The National Aeronautics and 
Space Administration (NASA), the Federal Aviation 
Administration (FAA), airport operators, and the airline 
industry are all interested in methods to improve airport 
capacity. The NASA is conducting a Terminal Area 
Productivity (TAP) Program to provide the necessary 
research to support the FAA and industry in safely 
achieving clear-weather (visual flight rules) airport 
capacity in instrument meteorological conditions (IMC). 
The TAP Program consists of four elements: Air Traffic 
Management, Aircraft-Air Traffic Control Systems 
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Integration, Low-Visibility Landing and Surface 
Operations, and Reduced Spacing Operations. 

The wake vortex research described in this paper is 
a part of the Reduced Spacing Operations element of 
the TAP Program. Improvements in landing frequency 
through reduced in-trail spacing have the potential to 
improve system capacity by KD-15%. 1 This potential 
capacity improvement has renewed interest in wake 
vortex research. 2 The spacing required to avoid the 
wake turbulence of the preceding airplane is one of the 
limiting factors in safely reducing in-trail spacing. 

The need for sufficient spacing between aircraft to 
insure an acceptable vortex encounter upset is not 
exceeded is most critical for aircraft near the ground 
during landing and takeoff. 3 The degree of upset 
mainly depends on the relative size of the vortex 
generating and vortex penetrating airplane, and the 
extent the wake decay imparted by the atmosphere. The 
vortex initial energy or strength and the resulting 
response to a wake are directly related to aircraft size 
and weight. 

The knowledge base on which current wake-vortex- 
imposed separation standards are based evolved during 
the late 1960's and the 1970’s. Currently, the U.S., 
aircraft are classified as "Heavy" (300,000 lb. or 
greater), "Large" (between 12,500 and 300,000 lb.), and 
"Small " (less than 12,500 lb.) based on maximum gross 
takeoff weight. 3 ^ 4 Although technically a large aircraft 
by weight, the B757 is currently given a special 
classification when it is the lead aircraft. The 
separation distances are based only on the 
classifications of the leading airplane and the following 
airplane. No account is made for factors such as 
airplane configuration or actual weight, wing span, 
atmospheric turbulence and winds, or local terrain 
effects. Therefore, the ability to accurately model the 
wake hazard and determine safe separation distances for 
a wide range of aircraft, weather, and operational 
scenarios may provide the basis for significant increases 
in airport capacity. 

One of two conditions must be met to safely reduce 
the spacing between aircraft. The wake of the 
preceding aircraft must have either transported out of 
the intended flight path of the following airplane or it 
must have sufficiently decayed such that the 

encountering airplane can safely land. It is the 

determination of what constitutes sufficient decay that 
this study will address. This simulation study tested 
various safe landing criteria to define the wake decay 
required for a B737-100 airplane to safely encounter a 
wake and land under the guidance of an auto-land 


system. The wake strength and encounter geometry 
were varied to establish a safe encounter boundary. A 
sensitivity study was also conducted to assess the effects 
of modeling methods and accuracy. 

Previous Wake Encounter Simulation Research 

A number of vortex-encounter simulations were 
performed during the 1970s and 1980s. Among these, 
three studies stand out as particularly relevant to this 
simulation. In 1974, Nelson and McCormick 6 studied 
vortex encounters with a batch simulation using 
analytical transfer functions to represent the pilot. 
Sammonds, Stinnet, and Larson 8 used a piloted wake 
encounter simulation to establish hazard criteria from 
pilot opinion. Hastings and Keyser 9 used a piloted 
simulation to study the effect of vortex decay on the 
initial response of a twin-engine transport. Each of 
these simulation studies used strip theory to model the 
wake effect on the airplane. The same method was used 
in this study and will be discussed later. These 
simulations differed from this study in that they were 
primarily focused on defining what constitutes a 
hazardous wake encounter. The goal of this study is to 
assess the sensitivity of candidate criteria to define an 
acceptable or safe wake encounter. In particular, the 
encounter must be weak enough that the airplane can 
continue the approach and landing without undue upsets 
to the passengers and crew. 

B737-1Q0 Wake Encounter Simulation 

The baseline B737-100 simulation was a batch 
version of the six-degree-of-freedom real-time 
simulation of NASA Langley’s Advanced Transport 
Operating System research airplane. The baseline 
simulation was modified to include the wake model and 
the strip-theory calculation of the vortex-induced forces 
and moments. Instead of using a pilot model, the 
control inputs were generated from an auto-land 
system. 5 The maximum control authority available to 
the auto-land system was less than the full control 
authority of the airplane. This is typical of most 
autopilot systems. The maximum roll-control available 
to the auto-land system was 37.4% of the full roll- 
control authority of the airplane. This may result in 
acceptable encounter boundaries that are considered 
overly conservative to a pilot flying the landing. 
However, if automatic landings are to be performed 
under reduced separation conditions then acceptable 
encounter criteria must reflect the limits of such a 
system. 
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Figure 1. Wake encounter simulation geometry. 



Figure 2. Two-dimensional wake model. 


The geometry of the encounter simulation is shown 
in figure 1. The airplane was initially trimmed for a 3° 
glideslope approach to landing at 137 knots and flaps 
set at 40°. It began 24,000 feet out from the runway 
threshold, at an altitude of 1,323 feet, and was correctly 
centered on the localizer and glideslope. A pair of 
vortices stretched from 24,300 feet before the runway 
threshold to 10,000 feet after the threshold. The 
altitude of the wake ( h w ) and the lateral placement (y w ) 
was fixed for each encounter simulation. The wake 
altitude was varied from 1200 to 25 feet. The lateral 
placement of the wake was also varied to determine the 
location at which the landing criterion was most 
sensitive. 


Wake Model 

An empirically derived two-dimensional wake 
vortex model was used to describe the wake of the 
generating airplane. This wake model was proposed by 
Burnham in reference 2 to fit field measured data of 
airplane wakes. The model was defined in the inertial 
axis system by the circulation (T), core radius (r c ) and 
location (y> z) of two counter-rotating vortices, as shown 
in figure 2. 


The wake model defines the tangential velocity ( Vq ) 
of a single vortex as: 
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where r is the radius from the center of the vortex. 
The sidewash (v) and downwash (w) velocity 
components at a point P are obtained by summing the 
contributions of the left and right vortices. 


v 


1 

r 

\ 

z p ~ Z R 


\ 

Z p~ Z L 


2n 

1 R 

r 2 +r 2 

r 2 +r 2 



- 

C R R ) 


V C L L J 

- 


( 2 ) 


- 

\ 


v 

r* 

y p ~ y l 

-r* 

y P -yR 

F L 

r 2 + r t 2 
K C L L J 

r 2 +r 2 

[ r <R +r * )\ 


( 3 ) 


It should be noted that there are some subtle yet 
significant differences between the velocity profile of 
the vortex model used in this study and the velocity 
profile of a Rankine vortex model. The tangential 
velocity of a Rankine vortex is defined as: 


Wr/lnr 2 ,r < r c 
1 Tjlnr ,r > r c 


( 4 ) 


The difference in the vortex velocity profile of the two 
models is illustrated in figure 3 for the same core radius 
and circulation strength. A Rankine vortex has twice 
the velocity at the core radius as the model proposed by 
Burnham. At large radii the models converge. The 
Burnham model was selected for this simulation 
because it provided a better fit to the measured wind 
tunnel and flight data that the simulation will be 
validated against. The effect of different vortex models 
on the simulation results is not part of this study but is 
part of future planned research. 



Figure 3. Tangential velocity profiles of two different 
vortex models. 
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For this simulation study the axis system of the wake 
model was aligned with the runway axis system shown 
in figure 1 . The lateral spacing and core radius of the 
vortices were set to represent the wake similar to a 
B727 size airplane. The core radius for both vortices 
was set at 2 feet (r = r =2). The circulation, 

altitude and lateral displacement of the vortices were 
varied throughout the study. The vortex pair were 
constrained to be symmetrical (zl = Zr . > Tl = T*). The 
lateral displacement of the wake relative to the runway 
centerline was referenced from the left or port vortex 
(y„ = yi ). The right or starboard vortex was space 84 
feet to the right of the left vortex (y R = y L + 84). The 
effect of the ground plane at the lower altitudes was not 
accounted for in the wake model. 

Vortex Encounter Aerodynamics 

The aerodynamic effect of the wake on the 
encountering airplane was modeled using a strip theory 
method similar to that used in references 6, 7, 8 and 9. 
Strip theory is a simple method in which the wing, the 
horizontal stabilizer and the vertical stabilizer are 
divided into a series of chordwise strips, as shown in 
figure 4. Each strip is treated as a 2-dimensional (2D) 
airfoil for which the lift at the quarter-chord point is 
computed as a function of the flow incidence angle at 
the three-quarter-chord point. The incremental 
contribution of each strip is summed to determine the 
forces and moments on the airplane. 

Each strip is defined by its area, 2D lift-curve slope, 
angle of incidence, dihedral angle, and the body axis 
coordinates of the quarter and three-quarter-chord 
points at the mid-span of the strip. The 2D lift-curve 
slope of each strip on the wing and horizontal stabilizer 



was weighted to yield the proper span load distribution 
and total lift curve slope of the airplane in a uniform 
flow field. The lift-curve slopes of the vertical 
stabilizer strips were weighted to yield the correct value 
of C Y p for the airplane. This resulted in an over 

estimate of and C T ^ . The significance of this over 

estimate will be discussed in the modeling error 
sensitivity section. A more detailed description of the 
strip theory force and moment calculation is provided in 
reference 1 1 . 


The change in the aerodynamic forces and moments 
due to the wake was determined by taking the difference 
between the strip theory calculation with the wake 
velocities included and the same calculation with the 
wake velocities set to zero. For example: 


w ^^wake wake 


( 5 ) 


The change in the force and moment coefficients 
due to the wake were then added to the baseline six- 
degree-of- freedom simulation coefficients. Because the 
wake effect is computed as the difference between the 
two strip theory calculations it is more important to 
correctly model the “slope” terms (e.g., C^and C y ) 

than the total force coefficients (e.g., C L or C y ). 


Landing Criteria 


As was discussed in the introduction, the objective 
of this study was to investigate the sensitivity of various 
landing criteria in defining the wake decay required for 
a B737-100 to safely encounter a wake and land. This 
analysis requires the establishment of some criteria to 
define a safe and acceptable landing. Several metrics 
were selected as candidate safe landing criteria based on 
the wake hazard research conducted in the 1970’s. 7, 8> 9 
Maximum bank angle and lateral and vertical deviation 
from the flight path were selected as performance 
criteria. The maximum bank angle limit was varied 
from 7° to 12° in 1° increments. 

A proposed auto-land certification requirement for 
satellite navigation systems (ref. 10) was used to define 
the lateral and vertical deviation limits. These limits are 
presented in Table 1 as a function of altitude. The 
dimensions correspond to approximately a 1° localizer 
deviation and a 0.3° glideslope deviation. The 
horizontal deviations were computed relative to the 
localizer centerline. The vertical path deviations were 
computed relative to a baseline vertical profile. The 
profile was basically a 3° glideslope that included the 
flair maneuver prior to touchdown. 
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Altitude 

feet 

| Lateral Limit 
feet 

Vertical Limit 
feet 

0 

! ±27 

0 

50 

±51 


100 

±75 

±15 

200 

±110 

±32 

250 

±118 

±36 

300 

±125 

±40 

400 

±158 

±51 

500 

±192 

±62 

750 

±275 

±89 

1000 

±358 

±116 

1250 

±442 

±143 

1500 

j ±525 

±170 


Table 1 - Proposed flight path deviation limits for a 
satellite based auto-land system. 


The only non-performance criteria examined in this 
study was the ratio of the wake induced rolling moment 
to the maximum roll control power of the airplane. 

‘max wheel 

This criterion was suggested in reference 12 from 
results of wake encounter flight tests. The maximum 


roll-control power ratio limit (8 w ) was varied from 0.5 
to 0.7 in 0.1 increments. It should be noted that the 
maximum roll-control power ratio for the auto-land 
system is: 




0.374 C 


■ = 2.6745^ (7) 


^max wheel 


Simulation Analysis 

A series of wake encounter simulations were 
conducted to derive the maximum wake circulation 
strength (T) for a given wake position (y WJ z w ) that 
satisfied each safe landing criterion investigated. For 
example, figure 5 shows the maximum circulation 
contours for a 10° maximum bank angle limit as a 
function of wake altitude and lateral position. The wake 
position where the bank limit is reached with the lowest 
circulation strength occurs when the left vortex is 
between 10 and 16 feet to the right of the airplane 
centerline and the wake is no more than 200 feet above 
the runway. Similar contours were created for each of 
the criteria investigated. The minimum wake strength at 
each altitude is obtained from the contours to define the 
maximum circulation boundary as a function of altitude. 
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Figure 5. Maximum circulation contours relative to the left vortex position for a 10° maximum bank angle limit. 
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Figure 6. Boundary of wake circulation as a function of wake altitude for various acceptable landing criteria. 

Figure 6 shows the wake strength boundaries as a 
function of wake altitude for the various acceptable 
landing criteria. The control power boundaries are 
presented as labeled dashed lines. The value in 
parentheses is the autopilot control power ratio. The 
vertical tunnel limit was not a limiting criterion relative 
to the horizontal limit and is not included in the figure. 



12° Limit 
11° Limit! 


At the lower altitudes ( h w < 150 ft) the horizontal tunnel 
limit is the limiting criteria. For wakes above 150 feet 
the control power ratio or maximum bank angle is the 
limiting criteria. The control power boundaries and the 
maximum bank angle boundaries are very similar. The 
control power boundaries are nearly constant with 
altitude, as are the bank angle boundaries between 200 
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Figure 7. Wake location range of minimum circulation for various acceptable landing criteria. 
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Figure 8. Sensitivity of modeling error and methods on 10° bank angle boundary. 


-o- Baseline 1 0° bank limit 
o 90 % of modeled wake effect 
D-* 110% of modeled wake effect 
-a- Wake effect on wing & horz. stabilizer only 
-v- Wake effect on wing only 



feet and 1000 feet. The large change in the bank angle 
boundaries below 200 feet is primarily due to the 
increase in the tracking gains as the localizer and 
glideslope signals narrow near the threshold. At 
altitudes above 1000 feet the variation may be due to 
reduced tracking gains or due to some of the wake 
effect being compensated for when the airplane is 
initially trimmed prior to the start of the run. The large 
shift between the 8° and 9° boundaries occurs when the 
maximum bank angle of the wake induced upset is less 
than the maximum bank angle of the recovery 
maneuver. 

Figure 7 shows lateral range of wake positions 
(position of left vortex) where the strength was a 
minimum for each criterion. This minimum would be 
considered the maximum acceptable circulation for that 
criterion. Although the control power boundaries and 
the maximum bank angle boundaries are similar 
functions of altitude, the lateral wake locations at which 
the boundaries are defined are slightly different. The 
wake location at which the control power and horizontal 
tunnel limits are reached with the minimum circulation 
occurs with the wake closer to the centerline of the 
airplane than the bank angle limits. The wake location 
of minimum circulation also moves farther outboard as 
the bank angle limit is increased. However, the 
difference in the location of minimum circulation is well 
within the variation in tracking performance of airplanes 


flying the localizer and is therefore probably not 
significant. 

Modeling error sensitivity 

A sensitivity study was conducted to assess the 
effects of modeling methods and accuracy. The 
analysis was conducted in two parts. The first part 
assessed the effect of mismodeling the wake encounter 
aerodynamics. This was done by altering the output of 
the wake encounter model by ±10 percent. The 
modified simulation was then used to establish new 
limit boundaries. Figure 8 shows the results of the 
±10% modeling error on the 10° bank angle boundary. 
The modeling error resulted in a 10% shift in the 
maximum circulation. This makes sense, in that the 
wake velocities are proportional to the circulation and 
the change in the aerodynamics is proportional to the 
wake velocities. A 10% change in one should result in 
a 10% change in the other. The effect of the modeling 
error on the lateral deviation boundary yielded the same 
results. 

A number of wake encounter studies have been 
conducted which only accounted for the wake effect on 
the wing and ignored the tail effect. The second part of 
the analysis assessed the sensitivity of this modeling 
method. Figure 8 shows the results of modeling just the 
wing and horizontal stabilizer, and the wing alone, on 
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Figure 9. Test of wake effect on a B737-100 airplane in 
NASA Langley’s 30x60 tunnel. 


the 10° bank angle boundary. The results indicate that 
the vertical stabilizer can have a profound effect on the 
boundary values. However, the vertical stabilizer effect 
may be exaggerated due to the current strip theory 
implementation. The current implementation assumes 
that the vertical stabilizer is the sole source of C Y ^ for 

the airplane and neglects the fuselage contribution. 
This results in the strip theory over predicting the 
vertical stabilizer contribution and yielding high values 
for C ; and C T ^ . A series of flight and wind tunnel 

tests to measure the wake effect on the airplane 
aerodynamics have been conducted at NASA’s Langley 
Research Center (figure 9). The results of these tests 
are currently being analyzed and should help resolve 
this issue. This may change the magnitude of the 
circulation values for some of the boundaries, but the 
general trends are not expected to change. 

Concluding Remarks 

This study is part of a planned series of studies to 
establish criteria for acceptable wake encounters. The 
control input for the simulation was generated from an 
auto-land system that had less than half of the full roll- 
control authority of the airplane. This may result in 
acceptable encounter boundaries that are considered 
overly conservative to a pilot flying the landing. 
However, if automatic landings are to be performed 
under reduced separation conditions then the acceptable 
encounter criteria must reflect the limits of such a 
system. 

The major results of this simulation study can be 
summarized as follows: 


a) The vertical deviation boundary was not a 
limiting criterion and the horizontal boundary 
was only limiting for wakes at 150 feet or below. 

b) The control power boundaries and the maximum 
bank angle boundaries are very similar. 

c) The lateral position of the wake where the 
acceptable landing criterion is exceeded with the 
lowest circulation occurs when the left vortex is 
offset to the right of the airplane centerline or 
conversely when the right vortex is offset to the 
left of the airplane centerline. 

d) The wake lateral position of minimum circulation 
is farther outboard for the maximum bank angle 
boundaries than for the control power or 
horizontal deviation limits. 

e) Sensitivity study results indicate that the vertical 
stabilizer can have a profound effect on the 
boundary values and must be accurately modeled 
to establish valid boundaries. 

Research is ongoing to validate the simulation 
encounter models with flight test and wind tunnel data. 
This may change the magnitude of the circulation values 
for some of the boundaries, but the general trends are 
not expected to change. 
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